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Abstract 

qj I We are interested in the structure of the Lcc vertex in the Yang-Mills theory, where c 

,^ . is the ghost field and L the corresponding BRST auxiliary field. This vertex can give us 

information on other vertices, and the possible conformal structure of the theory should 

- be reflected in the structure of this vertex. There are five two-loop contributions to the 

rN I Lcc vertex in the Yang-Mills theory. We present here calculation of the first of the five 

^ I contributions. The calculation has been performed in the position space. One main feature 

of the result is that it does not depend on any scale, ultraviolet or infrared. The result is 

expressed in terms of logarithms and Davydychev integral J(l, 1, 1) that are functions of the 

ratios of the intervals between points of effective fields in the position space. To perform the 

calculation we apply Gegenbauer polynomial technique and uniqueness method. 
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1 Introduction 

Recently it has been shown that the effective action of the J\f = 4 SYM written in terms of the 
dressed mean fields does not depend on any scale, ultraviolet or infrared [1, 2]. The theory in 
terms of these variables is invariant conformally. Therefore, investigation of this action should 
be greatly simplified. It might be possible to fix all three-point correlation functions up to some 
coefficient. For example, the three-point function of dressed gluons in the Landau gauge could 
be found in such a way. Conformal symmetry does not help a lot in the case of the four-point 
function since an arbitrariness arises. However, in any case, it is important to check all these 
statements directly by the precise calculations of the vertices in terms of the dressed mean fields. 
Lcc vertex is the most simple object for this calculation, especially in the Landau gauge. In that 
gauge it is simply totally finite. This fact has been indicated first in Refs. [1, 3]. 

The purpose of the study is to calculate the Lcc vertex at two-loop level for AA = 4 super- 
Yang-Mills theory, which should confirm the statement of Refs. [1, 2] that were derived from 
the results of Refs. [4, 5, 6, 7, 8]. The present paper contains the calculation of the first needed 
diagram, and also all needed formulas. The evaluation of other diagrams is a subject of our 
future investigation. In our opinion, it is the first calculation of two-loop three-point diagrams 
in general kinematics for a real physical model. The calculation is very nontrivial and can be 
used by others in some different studies. 

With a little modification, the investigation can be applied also for the finite (and of course, 
singular part) of the vertex at two-loop level for nonsupersymmetric theory. As to the singular 
part in higher orders for nonsupersymmetric case, in Ref. [9] the explicit three-loop computation 
of the anomalous dimension of the operator cc has been carried out. The vertex Lcc is convergent 
superficially in Landau gauge in any gauge theory, supersymmetric or nonsupersymmetric. This 
fact is a consequence of the possibility to integrate two derivatives by parts and put them 
outside the diagram on the external ghost legs due to transversality of the gauge propagator. 
If we make a change of the normalization L to Lg, as it takes place in Refs. [9, 10], than the 

superficial divergence would mean that Z^ZgZc = 1. (Renormalization constant of 5 is Zg, and 

1/2 
the renormalization constant of c is Zc ). This is the condition of superficial convergence of this 

vertex. Formally, this result holds to all orders of the perturbation theory due to the so-called 

Landau ghost equation, stemming from the fact that in the Landau gauge the Yang-Mills action 

is left invariant by a constant shift of the Faddeev-Popov ghost c (Ref. [11]). Note that it is 

not like the normalization used in Refs. [1, 2] where Z^Zc = 1. Moreover, the renormalization 

of ghost and antighost fields in our Refs. [1, 2] was chosen to be independent. This is different 

from Ref. [9] in which Zc = Z^- Thus, in terms of conventions of Ref. [1, 2] Z^ = Zc = 1. 

It means the field c does not get the renormalization in our convention. This coincides with 

the old results obtained from the antighost equation of Ref. [11]. Neither does the external 

field L get the renormalization in this convention. Starting with the two loop order infinities 

reproduce the renormalization of the gauge coupling in nonsupersymmetric theory. If we change 

the convention of the renormalization, the relation between the renormalization constants of 

paper [9] for nonsupersymmetric case can be reproduced. 

Knowing the structure of the Lcc vertex, one can obtain other vertices in terms of this one 

by using Slavnov- Taylor identity [12, 13, 14, 15, 16, 17] which is a consequence of the BRST 

symmetry [18, 19]. Moreover, the algorithm for obtaining these structures is expected to be sim- 



pie, due to the simple structure of the Lcc vertex, in particular due to its scale independence. 
Similar arguments can be applied to A/" = 8 supergravity [20, 21], and to other theories which 
possess high level supersymmetry to garantee good properties of the correlators. In some theo- 
ries, for example Chern-Simons field theory, nonrenormalization of gauge coupling is protected 
by topological reasons and similar approach is valid near the fixed points [22] in the coupling 
space. 

The article has the following structure. Section 2 contains the derivation of gluon and 
ghost propagators in the position space. A review of the one-loop results and all two-loop 
diagrams, contributing to the problem, are given in Section 3. Section 4 demonstrates a useful 
representation for the two-loop diagram (a) which is the subject of this study. In Section 5 
we show basic formulas for calculation of the considered Feynman integrals. The calculation of 
diagram (a) is performed in Sections 6 and 7. Moreover, details of the calculations can be found 
in Appendices A and D. The most complicated Feynman integrals are evaluated in Appendices 
B and C. Section 8 contains conclusions and a summary of the results, and discussion about the 
future steps. 

2 Landau and ghost propagators in the position space 

In the momentum space the gluon propagator in Landau gauge is: 
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where the case a = 1 corresponds to the free propagator in four dimensions. We will assume 
that the Wick rotation has been performed and will thus work in the Euclidean metric. The 
number of dimensions is D = 4 — 2e. We formulate the rules in the momentum space and then 
we go to the position space. The Fourier transform of the (1) to the position space can be done 
with the help of the following formulas [25] 
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Thus, the propagator is of the form 
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with b = D/2 - a. 

The transversality condition 



a„(^-c- 



allows to determine the coefficient c 
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Thus, the free propagator in the Landau gauge is: 
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The ghost propagator in the momentum space is 
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In the position space the Fourier transform is 

ajd^ple- = 2^-vr^/^a(l)5,^ = 2^-\^l\{l){e - 1) 



2x, 



„2\2-e- 



(2) 



3 Diagram contributions 

The one- loop contribution to the Lcc correlator corresponds to the diagram of Fig. 1. This is 




Figure 1: One- loop contribution to the Lcc vertex. The wavy line represents the gluon propagator, the straight 
hnes are for the ghosts. 

the only possible one-loop contribution, since the ghost field interacts only with the gauge field. 
This one-loop result is simple since it does not require the integration in the position space and 
is proportional to the following expression: 
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(Xi - X2)'^(xi - X3)4 (Xi - X3)4(X2 - X3)4 (xi - X2)'^(x2 - X3)'^. 

As one can see, the one-loop contribution is simple, but it does not have a structure that 
is expected from conformal field theories ^ [23, 24]. This is because the external field L does 
not propagate, it is not in the measure of path integral. However, this vertex by Slavnov- Taylor 
identity can be related to the three-point function of dressed mean gluons and they are expected 
to have simple structure at least for the connected function in the Landau gauge in A/" = 4 
super- Yang-Mills theory. For this reason, it is important to calculate the next order of Lcc 
vertex since poles disappear there. Poles do not disappear in the triple correlator of gluons since 
they must be absorbed into the dressing functions of the gluons. 

Two-loop planar correction to Lcc vertex can be represented as combination of five diagrams 
depicted in Fig. 2. 
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Figure 2: The two- loop diagrams for the Lcc vertex. The wavy hnes represent the gluons, the straight lines the 
ghosts. Black disk in diagram (d) stands for the total one-loop correction to the gluon propagator. 



4 Integral structure 

As it was noted in the Introduction, in the present paper we analyse diagram (a) only. The 
derivatives on the ghost propagators at the points X2 and X3 can be integrated outside the 
diagram and be put on the external legs. Indeed, the result for the diagram (a) can be represented 
as the derivative 



J_/ (23^23), 

[23] V^'' ^ [23] 
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to an integral T which contains numerators of the propagators independent of X2 and X3. We 
have introduced for the brevity the notation 

[yz] = {y- zf, [yi] = (y-a;l)^.... 



^I.K. thanks A. Jevicki for clarifying this point 



and so on. 

These derivatives will simplify the Lorentz structure of the wavy gluon line and reduce the 
integrand to the scalar structure. The corresponding contribution is of the following form: 



Note that it is the one-loop contribution multiplied by the last two denominators. 

The Lorentz structure can be simplified. Indeed, because the scalar product 2(x2;)^(yz)^ 
[xy\ + [yz] — [xz], we have 



1 [yz][yl] + [yz\[zl] - (2 - e)[yzf + (1 - e){[yl] - [zl]f 

2 [yzf-'[ylf-'[zlf-' 

Thus, the integral we have to calculate is 
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r= [ n n Hbl] + [vAi^^ - (2 - ^)[vz? + (1 - e){[ylf + [zlf - 2[yl][zl]) 

y ^ [yz]2-[yl]2-[zl]2-[y2]i-^[z3]i-^ ' ^> 

where we use the notation Dy = tt~ ''^d y. ^ 

We expect that the diagram T is finite in the limit e — > 0. Infrared divergences in the position 
space can be analysed in the same manner as it has been done for the ultraviolet divergences 
in the momentum space in the Bogollubov-Parasiuk-Hepp-Zimmermann /^-operation. From the 
expression above, for example, it can be seen that the infrared limit |x| ^ oo is safe in the 
position space in all the subgraphs and in the whole diagram. In the ultraviolet region of the 
position space each of the integrations is safe, too. 

Since the diagram is finite, does not matter where precisely the "e" is. In a certain sense, 
it is possible to change the indices in the propagators by adding multiples of e. In this way 
we can achieve the possibility to use the uniqueness relation [26] to calculate at least one of 
the two integrations by the bootstrap technique. Deviations in logarithms in the integrands, 
after changing the indices, cannot change the results in the limit e — > since they present finite 
construction times e. 

Thus, with the accuracy 0{€^), the integral above can be transformed as 

r, r, b^]bl] + [yz][zl] - (2 - e)[yz]^ + (1 - e){[y^? + [^1]' " 2[|/1][^1]) . . 

y H2[yl]2-2e[^l]2-2e[y2][z3] ' ^^^ 

where e = — 2e/(l — 2e), that corresponds to c = —2/(1 — 2e) and/or 6 = 1. This change e — > e 
is necessary, since we need to keep transversality in the position space after changing the index 



^we have introduced new D-dimensional measure Dx = Tx^~d^x. That choice helps us to avoid the nor- 
mahzation factor 7r^"/(27r)^ = l/(47r)^", coming in calculation of every loop, with the standard measure 
(f X / {2-k)^ .The above factor -k^'^ comes from integration itself. So, with the new measure Dx the final results 
in n-order of perturbation theory should be added by the factor Q"/(47r)^"", where as = gl is the coupling 
constant. Thus, in perturbation theory with the expansion parameter as/(167r^), the results in the new measure 
Dx do not obtain an additional factors in MS like scheme. 



of gluon propagator from 1 — e to 1, according to the formulas of Section 2. Transversality is 
necessary requirement to avoid problems with ultraviolet divergence in the position space. 

Moreover, the diagram is symmetric with replacement {y, 2} <-^ {z, 3} and, thus, we replace 
Eq. (6) by 

r=j[ n n M[^i] - (i -^/2)H^ + (i - 6)([zi]^ - [yi][zi]) \ 

|{y,2}^{z,3}| . (7) 

Following to the Eq. (3), the final results for the first diagram V in Fig. 1 can be represented 



as 



5 Technique of calculation 

To calculate expression (7), we need to use uniqueness method [26, 27, 25] and Gegenbauer 
polynomial technique (GPT) [28, 29]. Let us to give a short review of the uniqueness method. 
The GPT will be used only for the most complicated diagrams in Appendix A. All needed 
formulas for the GPT application can be found in [29]. 

The uniqueness method contains several rules to calculate massless chains and vertices alge- 
braically, i.e. without a direct calculation of D-space integrals. 

1. The results for chains J{ai, 02) have the form 
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J{ai,a2) = J Dxj-^^^^^^-^=A{ai,a2,as)j-^ , (9) 

where 

A{ai,a2,a3) = a{ai)a{a2)a{a^), Dx = 'K^~d x 

and 03 = D — a2 — cti and dj = D/2 — Oj. The point xi can be shifted to xi = 0. We have 

— — n 
introduced new D-dimensional measure Dx = vr 2 d x. Note that chain can be considered as 

the vertex with one propagator having power 0, i.e. J{ai,a2) = J(ai,a2,0). 

2. Uniqueness method [26, 27] (see also nice review [25]): if ai + a2 + as = D, then 

•^("1' «2, as) - / ^- [,!]«, [J].. [^3].3 = ^("1' «2, »^\uMi3M23f^ ' ^'"^ 

3. Integration by parts procedure (IBP) [26, 27, 25].^ 



'in the momentum space the equal relation for triangle is also very popular procedure (see [30]). 



Including in the integrand of J(ai, a2) 03) the function 3^(xl)^ and applying integration by 
parts, we obtain 
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The first term on the r.h.s. is equal to zero. Performing the derivative in the second term, after 
little algebra we obtain IBP relation 



(D — 2a\ — a2 — 03] </(ai,a2,a3) = 02 (</(«! — 1, 02 + 1,03) - [12] J(ai,a2 + 1,03) 

+a3( J(ai - 1,02,03 + 1) - [13] J(ai,a2, "3 + 1)1 • flU 



The relation has symmetry with respect to 02 ^^ 03- The index a\ has a special role, and 
the corresponding propagator having the power a\ will be called distinguish line. 

6 Cancellation of poles 

Now, using results of Appendices A and B, we obtain the final results for expression T. For this, 
it is convenient to consider the following combination: 

- (l-2e)[12]i-^[13]i-^ ^ 
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The last identity holds because 

r(l + ae) = exp 
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where 7 and C,{k) are Euler constant and Euler numbers, respectively. Thus, we obtain 

2(1 -e) 
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For the second term T2 we have 
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where 
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Thus, we obtain 
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The third term is very simple 

T3 = ([12] + [13])j(l,l,l), 
where for integral J(l, 1, 1) we can use Davydychev formula (see [31]): 
2 
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where 
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Note that the results (13) have a clear symmetry {2 ^-> 3}. 
The last term T4 has the form 
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Consider terms in the brackets. Ah 0(1) terms are cancelled: 
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At 0(e) level, all logarithms are also canceled. Indeed 
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So, the terms in brackets are 
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Combination of three terms 

(Ti + T2 + T4) /r(l + e) = 2 - L - 2[23] J(l, 1, 1) + 0(e) = Ti + T2 + T4 

is finite and, thus, 

T = 2-L+ ([12] + [13] - 2[23]) J(l, 1, 1) . 

Because eyl(l, 1, 2 — 2e) = 1 + o(e), we obtain 
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7 Final result 

To obtain the results (8) for the diagram (a) of Fig. 1, we apply the differentiation procedure 
(3) to ther.h.s. of Eq. (15). 

To write it in the more convenient for the differentiating form it is convenient to represent 
the expression in terms of the function J(l, 1, 1). 

It is better to represent the integral T as three terms: 

T = 2T(^) - t(2) + t(^) t(^) = ^ 

[12] [13]' 

^'^ - p|l^ ^' ^^'^ - P^ ([12] + [13] - 2[23]) J(l, 1, 1). 

To obtain the final results for the expression V representing the first diagram (a), we apply 
the projector 
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to the integral T, i.e. 
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2. To obtain V'^'^' , we represent it as the sum of three terms: 
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The most complicated second part of V^'^' has the form 
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and, so. 
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Thus, we have 
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where 

B2 = ([12] - [13])2 - 2([12] + [13]) [23] + [23]^ = B^ . 
3. The third term V^^' can be also represented as the sum of three terms: 

(3) (3) (3) 

where the terms Vi , V^ and V^ can be found in Appendix D. 

Collecting all these terms together, we obtain, after some algebra 
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4. Now the result for the expression V representing diagram (a) has the form 
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where 

Ai = 251-4^2 + 254 = 8[12][13]+ 4([13] + [12]) [23]- 4[23]2, 

A2 = -B1 + 2B5 = [12]2-3[13]2+4[12][13] + (5[13]-[12])[23]-2[23] 
^3 = 1^2,2^31, 

^4 = S7 = 4[23]3-6([13] + [12])[23]2+ ('([13] + [12])^-2[12][13])[23] 
+ (([13] + [12])'-6[12][13])([13] + [12] 

8 Conclusions 

In this paper we have shown, among other things, that the first of the five two- loop contributions 
in the correlator Lcc does not depend on any scale. The calculation has been performed in 
the position space and in the Euclidean metric. For this particular contribution it is a direct 
consequence of the transversality of the gluon propagator in the Landau gauge. The same is 
true for the two other vertex- type contributions. The AA = 4 supersymmetry does not play 
any role in the scale-independence of the first three contributions. It is important only for the 
cancellation of poles between the other two contributions of propagator-type corrections. The 
present study was necessary to investigate the precise structure of the two-loop contributions. 
By the ST identity the Lcc correlator can be transformed to the correlator of three dressed 
gluons. It is natural to expect that the conformal symmetry of the theory fixes the correlator of 
this triple gluon vertex completely up to some coefficient (that depends on the gauge coupling 
and number of colours). Since it is expected that the structure of the correlators of three dressed 
gluons is simple in the position space, the precise structure of the Lcc vertex in the position 
space is also simple. There are also other motivations in favour of the expected simple structure 
of the Lcc vertex. 

Indeed, the results for V and T contain the terms with different values of the transcendental- 
ity level. It was shown in [32] that some values oi N = A SYM variables have only terms with the 
same value of the transcendentality level at any order of perturbation theory (see [32, 33, 34]). 
It is possible that the same property is applicable to the Lcc vertex. If this is the case, then the 
final two-loop result for the Lcc vertex should contain [in the numerator of the r.h.s. of Eqs. 
(15) and (16)] only J(l, 1, 1) vertex and/or ("(2) Euler number. 

Finally, we would like to mention another obvious consequence of the considerations made 
here and in Refs. [1, 2]. All this can be applied to any gauge theory with only one coupling 
(gauge coupling) whose beta function is vanishing at every order. It means that the correlators 
of the dressed gauge bosons in that theory do not depend on any scale in the transversal gauge. 
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9 Appendix A 

The Appendix A is devoted to evaluate the expression for T. 

The integrand in brackets on the r.h.s. of (7) is a sum of several parts. The result, coming 
from term ~ [yz]'^ in numerator, is 

'' = /^^^" [,l]2-2.[,i]U[,2][,3] =-^'(^''-'^'^^ [12]Ml3]^- • 
This integral is very simple and requires only the use of Eq. (9). The term ~ [yl][zl] leads to 

""-"^ WEiPlHP«iI = •"''■'■'- '^'IBF''<' - 3.. 1 + ., 1) = 

^(1,2,1 -2e)A(2-3e,l + e,l) ^ 



[12]i-2^[13]i-2^[23]2^' 
where Eqs. (9) and (10) have been used simultaneously. Because of the relations 

A(l, 2, 1 - 2e) = -(1 - 2e)A(l, 1, 2 - 2e), A{2 - 3e, 1 + e, 1) = " ^^^(2 " 3e, e, 2), 

X 2(l-3e) ,, X .. ^ 

A(l-3e,l + e,2) = ^-^^(2-3e,e,2) (Al) 

we obtain 

/4 = ^(l,l,2-2e)A(2-3e,e,2) ^^^^,_,^^j^,„,^^3^, . 

Using the uniqueness relation (10), we have for the term ~ [yz][zl] 

/■ 1 1 

h = DyDz , ,, ,^ ^ , ,, ^ , ,, , =^(1,1,2 -2e)-^-n — J(2-3e,e, 1). 

J ^ [yz][yl]2-2^[zl]i-2^[y2][z3] ^ ' ' ^[12]i-^ ^ ' ' ^ 

The last integral 

' DyDz 



H2[yl]2-2e[^l]-2e[y2][^3] 



can be reduced by integration by parts to the basic integrals. It is better to consider first the 
integral /i and to apply IBP to the vertex with the center in point z and the distinguished line 
with the power 1 — 2e. The result has the form 



2e/i = h-h+ I DyDz 



[zl] - [31] 



M[yl]2-2e[^l]l-2.[y2][^3]^ 
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Evaluating the integrals on the r.h.s., we obtain 

/3 = ^(1,1,2 -2e)- 



+ 



[12]i- 
^(1, 2, 1 - 2e)A{l + e, 1, 2 - 3e) + ^(1, 1,2- 2e)A{e, 2,2- 3e) 



2eJ(2-3e,e,l) - J(l-3e,e,2) 
1 



[12]i-2^[13]i-2^[23]2^ 



Because of the relations (Al), the integral I3 has the form 

1 



/3M(l,l,2-2e) 



2eJ(2-3e,e, 1) - J(l-3e,e,2) 
1 



[12]i- 
+2^(e,2,2-3e 



[12]i-2^[13]i-2^[23]2^ • 



(A2) 



Both r.h.s. integrals J(2 — 3e, e, 1) and J(l — 3e, e,2) have singularities at e ^ and it is 
convenient to express J(l — 3e, e, 2) through J(2 — 3e, e, 1) and J(l — 3e, 1 + e, 1). The last integral 
is finite at e — > 0. 

Applying IBP to the integral J(l — 3e, e, 2) with the distinguished line having the power 2, 
we obtain 



-J(l-3e,e,2) 



J(l - 3e, 1 + e, 1) - [23] J(l - 3e, 1 + e, 2) 



+ (1 - 3e) 



J(2-3e,e,l) 



-[13]J(2-3e,e,2) 



eJ(l - 3e, 1 + e, 1) + (1 - 3e) J(2 - 3e, e, 1) 



- eA{l + e, 2, 1 - 3e) + (1 - 3e)A(e, 2,2- 3e) 
Because of relations (Al), integral J(l — 3e, e,2) obtains the form 



1 



[12]-^[13]i-2^[23]2^ • 



-J(l-3e,e,2) = eJ(l - 3e, 1 + e, 1) + (1 - 3e) J(2 - 3e, e, 1) 



l-3e 
l-2e 



A(e,2,2-3e) 



[12]-^[13]i-2^[23]2^ ■ 



Then, for the integral I3 we obtain 

h/Ail, 1, 2 - 2e) = p|j3^ (1 - e) J(2 - 3e, e, 1) + eJ(l - 3e, 1 + e, 1) 

1-e 1 

+ 1^:2^^'' ^' ^ " ^'^[12]l-2e[13]l-2e[23]2. 

Combining all the results, we obtain for expression T: 

T = / + /(2 ^3) , 
where / = /i — (1 — e/2)/2 + (1 — e)(-^3 ~ I4) has the following form: 

"^"^^ ^ "^ (2-3e)J(2-3e,e, l) + eJ(l-3e, l + e,l) 

1 



(A3) 



^(l,l,2-2e) [12] 



l-e 



+ 



' ^(.,2,2-3.)l'^l"l'=" 



l-2e 



[23] 



2e 



:i-e)A(l,l,2-2e) 



[12]i-[13] 



l-e 



(A4) 
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This formula presents the result for the expression T. It contains two integrals J(2 — 3e, e, 1) 
and J(l — 3e, 1, 1 + e) which cannot be calculated by rules from the previous section and will be 
calculated in Appendix B by using GPT. 

10 Appendix B 

In Appendix B we consider the integral J(2 — 3e, 1, e). 

1. There is a relation between two integrals J(2 — 3e, 1, e) and J(A, A, 2A) where A = 1 — e. 
The last integral is more convenient to calculate by GPT. 
First, we transform using the uniqueness relation 



J(2 - 36, 6, 1) = T^Jr^Ji^ - 3e, e, 1)[23]'-'^ = J 



Dx 



1 



[23] 



l-2e 



Dx- 



1 



Dy. 



[23]i-2.[a.i]2-3.[a.2]^[x3: 
1 1 



[23] 



l-2e 



xl]2-3^ J ^ [yx]i+^[y2]i-^[y3]2-2e A{l+e,l-e,2- 2e) 
i_2. A(2-3e,l + e,l) f „ 1 



[23]^ 



^(l + e,l-e,2-2e) 



Dy 



[yl]i-^[y2]i-^[y3] 



l-eL,Ql2-2e 



i_2, A{2 36,1 + 6,1) _ 1 _ 2 - 26). (Bl) 

^ ^ ^(1+6,1-6,2-26) ^ ' ' ^ ^ ^ 

2. Thus, we calculate the integral appearing above, by applying GPT (one can put ^3 = 
by shifting the arguments), following [29] 



/oo 
Dx J2 MnW 
n=0 
-. 00 

— - y: M„(A)xf 

^ y^) n=0 



^^ll^l2.■■^ln^^^l^J■2■■■^J■ 



/i2---Mn^Ml/^2---Mn 



(xl) 9{lx) 



lA+n 



[1] 



\+n 



1 1 



^(211 



1 



[2]3A+n-i(3A + n-l)(l-2A) 
0(12) 0(21) /[l]^"+^~ 



1 0(12) 



1 



[l]3A+n-i(^ + A) \^3A + n-l) 1-2AV[2], 

1 / 0(12) 0(21) /[1]^"+^^ 



+ 



+ 



[2]2A-i [i]A+n(2A-l)(n-A + l) [l]3^+«-i(n + A) \^2A-1 n - A + 1 V [2] 
We introduce notation 

M„(A)^^^^^^, e{i2)^e{xl-xl). 

Here, M„(A) is the coefficient at the traceless products appearing in the Gegenbauer polynomials 

Cn{xiX2) 



C^{xiX2) = Mn{\)^ 



Ulfl2 ■■■^n ^fJ'lfJ'2 ■■■^n 



OCrt 



(x?)"/2(a;i)"/2 
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Transforming the previous expression one obtains: 



1 1 

J(A, A, 2A) = Y^XiJT^ ^ ^^(^)^MiM..../^..mM..../^. 



X2 

ra=0 



0(21) ' ^ ^ 



[2]3A+n-i (3A + n - 1) 



11 1 \ ,, , / 1 1 11 1 

+ 12 



[1]2A-1 [2]"+A„_A + i; ' ^ V[l]3A+n-l3_X + „_l [2]2A-1 [l]A+n „ _ ^ + ^ 

^ ^ V21)i2i(A,A,2A)+0(12)ii2(A,A,2A) 



r(A) (1 - 2A) 



where ji2(A,A,2A) =i2i(A,A,2A)([l] ^ [2]^ 

Thus, we can consider below only the case 0(21). It is convenient to use Gegenbauer poly- 
nomial itself to reconstruct the results for J(A, A, 2A). 

Note that we can re-present the expression for ji2(A, A, 2A) as 

j2i(A, A,2A) = E Q'(^i^2) (^ [2p-3e2 + ^_3e " [iF^W^^TTlj' 

where ^ = \/[l]/[2]. We would like to note that GPT has been used before only for the 
propagator-type diagrams (or for the vertex ones with very specific kinematics - see the last 
entry of Refs. [28]), where the results could be represented as some numbers, i.e., for example, 
with ^ = 1 and C^{1) in above formula. Here, GPT is applied for the first time for the diagrams 
having two independent arguments. So, we need a technique for reconstruction of the final re- 
sults from the expansion in Gegenbauer polynomials. To obtain it, we represent the above series 
in terms of integrals 

i.2i(A,A,2A) = 



n=0 



[i]i-2^[2]i-2^^ ;S) ^ 

1 /•« duj {uj^-^^ - uj^'^ 



[l]i-|^[2]^-i^ ./o (l-2cos0Cc; + Cc;2)A 
New notation is introduced 

(2^12^2) = COS0 . 

The second integral over w in the above expression can be re-written as 



5 w^"^ f^ 2 

duj- = / dujuj'^ {I — 2 cos 9uj + uj Y 



(l-2cos6la;-Fa;2)i-^ Jq 

/ dLO^-{l-2cos0uj + uj^Y + I diOLO 
Jo e Jo 



1 2 cos — uj 

+ 



UJ 1 — 2 cos Quj -|- ( 



^ , A,. . . 9^. /"^ , . 2cos6'-a; 

(\-2zo?,Quj^J^Y~^ 
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e Jo 



« , Acos 6u;' -3u;^+' 



'l-2coseuj + uj'^) 



2M-e 



1 / [l][12]n '^^ /-g 4cosgc^^-3c^^+^ 

7 \W~) ^ Jo (l-2cos^t^ + w2)i-e • 

In such a way we have extracted the pole in e. The integral in the last line is not singular. Now 
the total expression for J21 is 

i2i(A,A,2A) = 
[12]-- / 1/[12]\=^- [12]- /•« -4cosg^- + 3^i+-+u;i-3- \ 

[l]l-2e[2]l-2e 1^ e { [2] J ^ [l]e/2 [2]./2 jg «'^ (1 - 2 COS ^W + UJ^)!-. ^^ ^^^^ 

To check the self-consistency of the result, we have to check that the theta- functions disappear 
and that we left with a Lorentz-invariant structure. For the poles this is obvious. Let us check 
this for the zeroth order in e. The integral can be decomposed as 

«, -4cos6'a;- + 3a;^+-+a;^"3- /"« , -AcosO + Auj 

duj ; 7TZ-1 = / duJ 7T + 

(l-2cos6'cj + w2)i-^ io 1-2COS 6*^ + ^2 

/ , „ fi , Inuj l< , (-4cos6' + 4tj)ln(l -2cos6la; + a;2)\ 

e -4 cos 61 / dw ^+/ dto- — —^ ^ \ + o{e) . 

\ Jo l-2coseuj + uJ^ Jo l-2coseuJ + uj^ J ^ ^ 

Since we have to calculate the finite part of the initial diagram, the higher orders in e do not 
contribute in the limit e ^ 0. 

The first and the third integral can be easily calculated 

/■« -4cos6' + 4u; f^ , 9^ . 9^ [12] 

/ duj 0=2 dln(l-2cos9uj + uj^) = 2\n(l-2cos9e + C) = 2ln^ , 

Jo l-2coseuj + uj^ Jo ^ ^ J [2] ' 

<, (-4cos6' + 46j)ln(l-2cos6la; + a;2) /"« „ 2. . 2x , 2 [12] 

dco- — ^— ^ -= dln^(l-2cos9uj + uj^) = ln^^ . 

l-2coseuj + uj^ Jo [2] 

Moreover, by the conformal substitution to -^ 1/uj it can be immediately verified that 

« , Ino; /■!/? Ino; 

duj 7T = / duj K, , 

1 — 2 cos Ouj + uj'^ Jo 1 — 2 cos duj + o;^ 

i.e. the integral is symmetric under the exchange [1] ^^ [2]. 

Expanding the r.h.s. of (B2) in powers of e, we obtain that the terms of zero order in e are 
absent 

-21nM + 21„I3 = 0, 

|2] |2] 

an that the logarithms of the first order in e are 

, 2 [12] , [12] ^ [12]2 ^ [12] , [12] 
- In^ V-r + In V-r In hnh = In V-r In V-r • 
[2] [2] [1][2] [2] [1] 
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Thus, we present expression (B2) in the form 

[12]- 



i2i(A,A,2A) 



- + e In-f-rln-r-p -4cos6' / duj- — ■ — ^ 

e \ [2] [1] Jo 1 — 2 cos Ouj + cu'^ 



[l]l-2.[2]l-2, 

The result (B3) has completely symmetric form under the exchange [1] ^^ [2] 

J2i(A,A,2A)=ji2(A,A,2A) 
and therefore 



(B3) 



J(A,A,2A) 



1 



1 



[12]- 



r(l-e)2e-l[l]i-2^[2]i-2^ 



A 121 12 ^ ^ /•? Inc^ 

+e In V-r In V-r - 4 cos 6* / duj ^ 

I [2] [1] io 1-2 cos Ouj + cj2 



(B4) 



All the theta functions disappeared and the result is Lorentz invariant. Now we restore X3 
component by shifting back the space-time coordinates xi and X2- Taking this into account, we 
obtain 



cos 6* 
and introducing the notation 

dio- 



(X1X2) 



1 [1] + [2] - [12] 

[l]l/2[2]l/2 2 [l]l/2[2]l/2 



Inw 



/o 1 — 2 cos Oio + a;2 
we can write Eq. (B4) in the form 



= /(a:i,X2,X3) = /i23 



J(A,A,2A) 



[12]- 



r(l-e)2e-l[13]i-^[23] 



l-e 



1 / [12] [12] J13] + [23] - [12] ^ 
"7 + ' r M M " [13]V2[23]i/2 ^12=^ 



(B5) 



3. Now we want to express the integral /123 through J(l, 1, 1). It is convenient to start 
with the integral J(A, A, 2A — 1), with the purpose to use the results of the previous subsection. 
Applying GPT we calculate the integral (one can put X3 = by shifting the arguments) 



„ 00 
J{X,X,2\-1) = J DxY^MnW 



„^lltJ.2■■■^J-nr^^^l^J■2■■■t^n 



'X 



n=0 

00 



Xf 



(xl) 9{lx) 



^ ' n=0 

1 



U+n 



^(211 



[1] 



\+n 



1 



1 



[x2]^ [x]2A-i 

1 



[2]3A+n-2(3^ + „_2)(2-2A) 

^(12) 0(21) /[l]^"+^~ 



1 9(12) 



[l]3A+n-2(„ + ^) \^3^ + „ _ 2 2 - 2A V [2] , 

1 ( 0{12) 0(21) /[1]^"+^^ 



+ 



[2]2A-2 [i]A+n (2A _ 2)(n - A + 2) [1] 



3A+n-2| 



n + A) \ 2A-2 n- A + 2 V[2] 
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Transforming the previous expression one obtains 



1 1 

r(A) (2 - 2A) 



E MnWxl 



fllfl2---fJ.n^fJ.lt^2---IJ. 



Xn 



n=0 



0(21) 



1 



1 



+ 0(12) 



1 



1 



J(A,A,2A-1) 
1 1 

[2]3A+n-2 (3A + n - 2) 

1 1 1 



[l]2A-2 [2]"+An_A + 2y ' ' V [l]3^+"-2 3A + n - 2 [2]2A-2 [i]A+n „ _ ^ + 2 

^0(21)j2i(A,A,2A - 1) + 0(12)ii2(A, A,2A - 1) 



r(A) (2 - 2A) 



where ji2(A,A,2A — 1) = j2i(A,A,2A — 1)([1] ^-> [2]j. Thus, we can consider below only the 
case 0(21). As in the previous subsection, it is convenient to use Gegenbauer polynomial itself 
to reconstruct the results for J(A, A, 2A) obtained above. Applying the GPT technique one can 
represent the above expression as 



J21 



(A,A,2A-l) = £a 



A/ 



1x1X2) 



1 



e 



1 



e 



„=o ,[2]i-3^1 + n-3e [l]-2^[2]i-^ n + 1 + e 

We now use the GTP formulas to present series in terms of integrals 



7 . ^n 



1 1 ri 



X1X2) 



n=0 



[2] 1-3^^ 



l-3e 



dujuj 



n~3e 



1 



i2i(A,A,2A-l) 
1 f< 



[l]-2.[2]i-^i+^7o 

f 00 



duLo' 



n+e 



1 ri 



-3e , .n+e 



duj- 



(uj-^^ - UJ^) 



2\X 



[l]|-|^[2]|-r ^ (l-2cos0c^ + c^2) 



1 



4e / duj- 



Inw 



+ o{e 



[l]5-f^[2]5-f^ [ Jo 1-2cos6lj + lj'^ 

-4e/i23 + o(e 



i r T 1 
2 



[1]2[2] 
The result has completely symmetric form under [1] ^^ [2]: 

j2i(A,A,2A-l)=ii2(A,A,2A-l) 
and 



J(A,A,2A-1) 



1 



[1]^[2]^ 



-2/123 + 0(1) 



Now we restore the xa-dependence and obtain the relation 

2 



J(l,l,l) 



[13] 1/2 [23] 1/2 



123 



(B6) 



(B7) 
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The following consequences can be derived from here: 
1 . 1 . 



-/ 



123 



[13] 1/2 [23] 1/2 ^^■^ [12]V2[23]i/2 
Thus, we can write Eq. (B5) in the form 



132 



'123 



[13] 



1/2 



[12]V2 



(132 



J(A,A,2A) 



[12]- 



7 + K^"|i|^"|i| + (f''] + t23]-[i2])j(i,i,i; 



r(l-e)2e-l[13]i-2^[23]i-2> 

11 Appendix C 

In Appendix C we consider the integrals J(l — 2e, 2, 2) and J(2e, 1, 2). 
1. The integral J(l — 2e, 2, 2) can be represented as 

J(l - 2e, 2, 2) = 1 a(3)9f J(l - 2e, 2, 1) = 1 ^(1 - 2e, 2, l)9f af ^^^^ 



4e 



4e 



Because 



[12]i-^[23] 

n(3)o(3) [13]^ ^ 

^ ^ [12]!-^ [23]!+^ 



l+e 



8e(l + e) 



[13]^ 



(23)^(13) 



tiy^'jjfi 



4e(l + e 



[12]i-^[23]2+^ [12]i-^[13]i-^[23]2+^ 
[13]^ , [12]^ 



+ 



[12]i-^[23]2+^ [13]i-^[23] 



l-e[OQl2+e 



+ 4e. 



4e^ 



1 



[12]i-^[13]i-^[23]2+^ 
1 



[12]i-^[13]i-^[23]2+^ 



we have 



[l + e] 



[13Y 



+ 



[12Y 



J(l-2e,2,2)M(l-2e,2,l) 
1 



[12]i-^[13]i-^[23]2+^ ■ 



.[12]i"^[23]2+^ [13]i^^[23]2+' 
2. The integral J(— 2e, 1, 2) can be rewritten, using uniqueness, as 

J(-2M,2) = [12]^-/i^.--— i 



[12] 



l+e 



1 



_^ , Dv 

yl(2 + e, 1 - e, 1 - 2e) J [x3]2 ' 



Dx- 



1 



-2^[x2][a;3]2 
1 



[12] 



l+e 



A{2,l-2e,l) 



Dy, 



[yl]l-.[y2]2+e[y.^]l-2e 
1 



A(2 + e, 1 - e, 1 - 2e) J " [yl]i-^[y2]2+^[y3]i-^ 
,1+, ^(2,1-26,1) 



[12]^ 



^(2 + e,l-e,l-2e) 



J(l-e,2 + e,l-e) . 



(CI) 
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The integral J(l — e, 2 + e, 1 — e) can be treated in a similar manner 

J(2 - 3e, e, 1) = [131' / Dx-, — t-, — ? — ttt-, — tt — ^-r = 
^ ' ' ^ ^ W [xl]i-^[x2]2+^[x3]i-^[13]^ 

[13]^ ,,, /. ., / ^2;-l- / Dy- ^ 



A(l,l,2-2e)y [x2]2+^y "^ [yl][y3][yx] 
^^3^,^(2 + .,2 -2e,-.)r 1 



A(l,l,2-2e) 7 ^[yl][y2]2[y3] 

, ,,yl(2 + e,2-2e, -e) , 

13 ' ^ J ' V^^ 1,2, 1) . 

^ ^ ^(l,l,2-2e) ^ ' ' ^ 

So, for the initial diagram J(— 2e, 1, 2) we have 

J(-2e,l,2) = [12]i+'[13]V(l,2,l), 

because 

^(2,l-2e,l) ^(2 + e,2-2e,-e) ^^ 
^(2 + e,l-e, l-2e) ^(l,l,2-2e) 

Evaluation of the integral J(l,2, 1) is very simple. Following Ref. [35], we apply IBP to 
J(l, 1, 1) with different distinguished lines ^: 

{D - 4) J(l, 1, 1) = 2J(0, 2, 1) - [12] J(l, 2, 1) - [13] J(l, 1, 2), (C2) 

{D - 4) J(l, 1, 1) = 2J(1, 0, 2) - [23] J(l, 1, 2) - [12] J(2, 1, 1), (C3) 

{D - 4) J(l, 1, 1) = 2J(1, 2, 0) - [23] J(l, 2, 1) - [13] J(2, 1, 1). (C4) 

Considering the combination 

[23]- < {B2) > -[13]- < (S3) > +[12]- < {B4) >, 

where symbols < {B2) >, < {B3) > and < (BA) > represent the results of Eqs. (C2), (C3) and 
(C4), respectively, we have 



{D - 4) (^[23] - [13] + [12] j J(l, 1, 1) = 2(^[23] J(0, 2, 1) - [13] J(l, 0, 2) + [12] J(l, 2, 0) 

-2[23][12]J(1,2,1) 

The expression on the l.h.s. is negligible when D ^ 4. Thus 

^(1, 2, 1) = pjpj ([23] J(0, 2, 1) - [13] J(l, 0, 2) + [12] J(l, 2, 0) 

and 

[12]113]' 



J(l,2,l)M(l,2,l-2e) 



[23] 



[12]-' - [13]-' + [23]- 



(C5) 



'in momentum space, similar analysis has been done in [36, 31]. 
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12 Appendix D 

In Appendix C we calculate the most complicated term V^^' contributed to the first diagram in 
Fig.l . 

1. The first term V^ has the form: 

yf V^(l, 1, 1) = ^^'^^f^^^'^ p|l^([12] + [13] - 2[23] 



The derivatives generate 

[12113] (f''^ + f''^-'f''^)=PFl3] 



^?^^^'^TTWHt([12] + [13] - 2[23]) = ([12][13]g^. - 2[12](23)^(13), 



+2[13](23)^(12), - 2[23](12)^(13), . 



After simple algebra, we have 



where 



^1 " [12]2[13]2[23]'^^^'^'^^' 



Bs = ([12] + [13])2 - [12] [13] - 2([12] + [13]) [23] + [23]^ 



2. The third term V3 is 

(3) [12] [13] [23]^ _ 2p /5(2)^(3) j(i 1 1^ 
^3 [12] + [13] - 2[23] " ^^"^^ ^^'^'^'^ '^'^ •^^^' ^' ^^ 

2 (^9^)9f J(l, 1, 1) + (23),(23).9f a(3) J(l, 1, 1) 

The first term in brackets generates the additional factor 

2(2x)^(3x)^ ^ [2x] + [3x] - [23] 
[2x][3x] [2x][3x] 

in the subintegral expression of J(l, 1, 1), where x is the variable of integration. The term has 
the form^ 

^a(2)9f J(l - 2e, 1, 1) = J(l - 2e, 1, 2) + J(l - 2e, 2, 1) - [23] J(l - 2e,2, 2) . 

(2) (3) 
The second term can be calculated similarly: (23)^(23),^5/j dl J{1, 1, 1) generates 

4(23)^(2x)^(23)^(3x)^ 
[2x][3x] 



■'Hereafter we change J(l, 1, 1) -^ J(l — 2e, 1, 1) to have regularization and uniqueness of the "star". [25] 
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in the subintegral expression of J(l, 1, 1). 

By analogy with above case we obtain, after some algebra 

(23)^(23),aj2)43) j(i _ 2e, i, i) = j(i _ 2e, 0, 2) + J(l - 2e, 2, 0) - 

2J(1 - 2e, 1, 1) - [23] V(l - 2e, 2, 2) 

,(3) [12] [13] [23]^ 



Then 



K' 



^ [12] + [13] - 2[23] 
I [23] J(l - 2e, 2, 1) + J(l - 2e, 2, 0) - J(l - e, 1, 1) - [23] V(l - 2e, 2, 2) 



+<^2^3 



V^^^ + \vP ,2 ^ ?> 



Taking the result for J(l — 2e, 2, 2) from Appendix B and using 

1 



J(l-2e,2,0)=^(l-2e,2,l)pj^ 



J(l-2e,2, 1) =y4(l-2e,2, 1) 



[13]^ 



[12]!-^ [23] 



l+e ' 



we have 



t7{3) 



y3^^^ + J(l,l,l) M(l-2e,2,l 



[13]^ 



+ 



[12]i-^[23]^ [12] 



l-e 



^1 + e 



[13]^ 



+ 



[12]^ 



[12]i-^[23]^ [13]i-^[23] 
1-e 1 



+ e 



[23] 



l-e 



[13]i-[12] 



l-e 



[12] 1- [13] 



l-e 



1 + e 



1 + ln 



[12] 
[23] 



+ €■ 



[23] 



[12] [13] 



Then 



^3(3) + [vi^\2^ 3} + 2J(1, 1, 1))/A(1 - 2e,2, 1) 
, [12] [13] [12] V [23] y [13] V [23] y 



Because eA{l — 2e, 2, 1) = — 1 + o(e), we have 



[12] [13] 



F3(3) + |y3(3),2w3} + 2J(l,l,l) 
2 ([12] + [13] - [23]) + [12] In III + [13] In ||| 



Thus, 



V, 



(3)_2([12] + [13]-2[23]) 



[12]2[13]2[23] 



2 ([12] + [13] - [23]) + [12] In III + [13] In ||| 

-2[12][13]J(1,1,1) 



(Dl) 
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3. The most complicated second part of V^'^' has the form 



y(3) _ p 



(MU 



^ V[12][13] 

(2) 

The derivative d}i generates 



([12] + [13] - 2[23])) # J(l, 1, i)| + |2 ^ 3 



3,(2) 



(mis] ([^'] + [^'] - '['^^O = PFM K^''^ " '^''^) ^''^'^ " 'f''^^''^ 



and after little algebra we have 

^^-^?^fTT^([12] + [13]-2[23]) 



V[12][13] 



[12]2[13][23]^ 



$i(12), + l>2(23), 



where 



^1 = [23] ([13] - 2[23]) , ^2 = [13] ([13] - [12]) - [23] (3[13] + 4[12]) + 2[23]^ 
It is convenient to represent (12)jy as (12),^ = (13),^ — (23),^. Then 



^""*" (ppi "^'1 + >''! - '!'=">) = PPPPF '*'<'=^'" + *^P=""1 ■ 



^J■'^^^l 



where 



$2 = [13] ([13] - [12]) -4([13] + [12])[23] +4[23]2. 
Thus, the considered term V^ has the form 

^2^^^ [12]2[13]2[23]2 = 2 j[13] [-^1(13)^ + $2(23)^]4=') J(l, 1, 1)|+|2 ^ 3 

= 2 



V'o^^^ + <! V^2 ^\ 2 ^ 3 



where 



V^^^ = [13] [$iVFf ) + ^sVFf ^ 



^(3) 



,(3) 



Consider first W2 = {23)udl J(l, 1, 1). The derivative generates the term 



2(23)^(3x)^ 
[3x] 

in the subintegral expression of J(l, 1, 1). 



[23] + [3x] - [2x] 



Thus, W^^^ is 



.(3) 



W2 -^(1,1,1) = -J(l-2e,0,2) + [23]J(l-2e, 1,2) . 
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Using Eq.(Dl) we obtain 



(W^f -J(l,l,l))M(l-2e,2,l) 



[12]^ 



^ = ^lnM 

[13]i-^[23]^ [13]!-^ [13] [23] 

Because eA{l — 2e, 2, 1) = — 1 + o(e), we have 



W. 



[13] [23] ^ ' ' ^ 



(3) (3) 

The operation W^ = (13)iy9t/ J(l, 1, 1) generates the term 

2(13)^(3a;)^ [13] + [3x] - [Ix] 



[3x] 
in the subintegral expression of J(l, 1, 1). 



[3x] 



Thus, Prf ^ is 



VFf ) - J(l, 1, 1) = - J(2e, 1, 2) + [13] J(l - 2e, 1, 2) . 



Taking the integral J(2e, 1, 2) from Appendix B and using Eq.(Dl), we obtain 

Vf)-J(l,l,l))M(l-2e,2,l) = 



[12]113]- 
[23] 



Thus, 
Then 
and 

where 



" 1 /I 1 1 Y 

[[23]^ M12]^ [13]^ ' [23]^ /J 


[231 [13] 


-i"-M-i-^<- 


,1). 



V, 



(3) 



V, 



(3) 



M"^^^M-"^^^M^t^^^("^^"^^^^^'^'^^ 



^^^^|i|+^^^^|i|+^^^(^'^'^^ 



[12]2[13]2[23]2 
-^2 + f [12] - [13]) f [13] + [12] - 2[23]) = f [12] - [13]) (2[13] + [12] 



+ 



2([12] +3[13])[23] -4[23]^ 



L,2 ^ 3 



[13] ($1 + $2) + [12]| ($1 + $2) , 2 ^ 3] = ([13] + [12]) [([13] - [12])' + 2[23] 



3 [13] + [12] +2[12][13] [23] . 
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